Higher order Shapiro steps in ac-driven Josephson junctions 
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We demonstrate that the well known phase-locking mechanism leading to Shapiro steps in ac- 
driven Josephson junctions is always accompanied by a higher order phase-locking mechanism similar 
to that of the parametrically driven pendulum. This effect, resulting in a 7r-periodic effective po- 
tential for the phase, manifests itself clearly in the parameter regions where the usual Shapiro steps 
are expected to vanish. 



oo 
(N 



> 

m 
o 
o 

&\ 
O 
00 
On 

O 

CZ3 



Oh- 



X 



The harmonically driven pendulum has been studied 
extensively over the past decades in many different con- 
texts. The driven pendulum has been one of the key sys- 
tems in nonlinear dynamics due to its simplicity and rich- 
ness in nonlinear phenomena such as phase-locking 0-^] 
and chaos Another reason for studying the driven 

pendulum is that the pendulum equation is the most 
widely used model for superconducting Josephson junc- 
tions . Particularly phase- locking of pendulum mo- 
tion to an ac-perturbation has been extensively studied 
in the literature due to the general interests in synchro- 
nization of oscillators and specifically due to the related 
Josephson junction, where phase-locking of the voltage 
response to an ac-current can occur in certain regions 
of the parameter space p). There are two distinctively 
different ways of driving a pendulum; direct and para- 
metric drive. The direct (torsional) ac-drive, which we 
will study in this paper, is the most relevant for techno- 
logically interesting systems such as the Josephson junc- 
tion, and it has been found to produce efficient harmonic 
|j],U and subharmonic [[[o| phase-locking of great interest 
for, e.g., the voltage standard |11|. The parametric drive 
(e.g., pivot oscillations) has mainly been studied from the 
point of view of nonlinear dynamics. The usual theoret- 
ical technique for predicting and analyzing the ac-driven 
pendulum is to assume that the driving frequency is much 
larger than the natural oscillations of the unperturbed 
pendulum. Within this framework, one can develop time 
separation analyses separating a fast linear response to 
the ac-drive from the overall slow nonlinear behavior of 
the system. The main difference between the effects of 
the two driving mechanisms is that the parametric drive 
usually leads to a 7r-periodic effective potential for the 
slow nonlinear behavior while the direct drive 

leads to the usual 27r-periodicity JjJB]. We will demon- 
strate in this paper that the direct ac-drive may lead to 
a 7r-periodic effective potential in certain regions of the 
parameter space and that this result has consequences 
for the well known Shapiro steps in the current-voltage 



(IV) characteristics of ac-driven Josephson junctions. 
We study the pendulum equation in the form, 

4> + a<p + sin (f> = r\ + e sin fit , (1) 

where <j> is the pendulum angle relative to vertical (down) 
and a is a normalized friction coefficient. Time is normal- 
ized to r = yfgjl, where g is the gravitational constant 
and / is the length from the pivot to the pendulum bob, 
•q is a dc-torque, normalized to mgl, m being the mass 
of the bob, and the normalized frequency, ft, and ampli- 
tude, e, define the direct ac-drive of the pendulum. In 
the context of Josephson junctions, <p is t ne phase differ- 
ence between the quantum mechanical wave functions of 
the superconductors defining the junction, a is the nor- 
malized dissipation coefficient due to transport of quasi- 
particles, r\ and e are currents normalized to the critical 
current of the junction, and time is normalized to the in- 
verse Josephson plasma frequency, r = y / hC/2eI c , where 
h is Planck's constant, C is the device capacitance, and 
I c is the critical current. Voltage across the Josephson 
device is given by V = (f>h/2e, uniquely relating Joseph- 
son voltage to pendulum speed || . 

Let us write the phase in the following form, 



9 + &+~(t) 



(2) 



where 3(i) is a function that oscillates with frequency 
ft, £ is a constant to be determined, and 9 is a phase, 
(6) = 0. Inserting Eq. (|) into Eq. (§) we obtain, 

9 + 3 + a0 + a£ + aS + sin (9 + £t + 5) = r? + e sin fit . 

(3) 

We will choose 3 so that, 3 + aS = e sin fit, and we then 
obtain, 



S(i) 



: sin (Qt + 7) 



(4) 



where 7 = tan 1 (a /ft) is a constant phase. Equation 
(@) then takes the form, 
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a9 + sin (6 + £i + H) = n - a£ 



(5) 



|?7 - «C| < |Gism(e + 5 Gl ) + G 2 sm(2e + 5G 2 



(10) 



Rewriting Eq. ([!]) in the form of Eq. (H) demonstrates 
how the direct ac drive can lead to parametric effects 
similar to those of the parametrically driven pendulum 
equation. 

We can now apply three consecutive transformations 
that will lead to an equation in variables where time- 
dependent terms are of order il~ 3 . We describe the gen- 
eral procedure as a generalization of the analysis pre- 
sented in Ref. [[l7| (which goes back to Poincare, see, 
e.g., Ref. fTs)] ) in the Appendix. 

To average the equation we rewrite it as a system of 
two ordinary differential equations (ODEs) (l^] 

9=p (6) 
p = rj — at; — ap — A sin 9 — B cos 9 , 

where A and B are given by, 

A = cos (ft + S) 
B = sin (ft + 5) . 

Carrying out the procedure described in the Appendix 
three times and neglecting terms 0(fl~ 3 ) or higher we 
obtain 

e = p (7) 

P = r\ - ai - aP - G x sin (6 + 8 Gl ) - G 2 sin (26 + S Ga ) , 
where 

-J fe (r) , if f/fi = k integer 



Gi = 







otherwise 



(8) 



G„ = J _ v Jn(r)J! „(r) 

2 1 f22( n _I)2 



if f/n = k = \ 

otherwise 



(9) 



and 



fiV^ 2 + a 2 ' 

t^Gj = tan 
8q 2 = tan 



-rCB) 
{A) 

. x 2({g}_ 1 {A}_ 1 ) 

({B} 2 ^) - ({A} 2 _ 1 ) 



Brackets, (■••), denote time average and {/} is given by, 
{/} = /— (/), where / is a periodic function. The mean- 
zero antiderivative, {/}_ 1 , is defined as, 

{/}-! = / {f}dt, ({/}_!) = 0. 

It is now clear that phase-locking of the pendulum mo- 
tion can exist for values of r\ given by 



For G2 = this leads directly to the well-known Bessel 
function expression for the Shapiro steps in ac-driven 
Josephson junctions jlj and this result is correct up to 
order 0(£1~ 2 ). However, we find that, for G\ = 0, 
i.e., at every node of the Bessel function we have 
phase-locking originating from the coefficient G2 to the 
7r-periodic term in Eq. (|7|). As a consequence, we may 
predict that the Shapiro steps in the IV characteristics of 
ac-driven Josephson junctions do not vanish for parame- 
ter values given by Jfc(r) = 0. 

In order to demonstrate this we have performed nu- 
merical simulations of the equation (|l|) and measured the 
ranges of phase-locking in n as a function of Q, e, k, and 
a. Figure la shows a typical normalized IV (rj, (</))) char- 
acteristic of an ac-driven system and we use this to define 
the magnitude of the locking range, Ar/k- The system pa- 
rameters have here been chosen to a = 0.3, f2 = 3, and 
r = 5. In Fig. lb we show the magnitude of the locking 
range, Aryj,, against e for ft = 3, a = 0.1, and k = 1. The 
solid line represents the usual Bessel function prediction, 
A771 = 2|Gi|, and the markers represent the numerical 
simulations. It is clear that as long as G\ (the relevant 
Bessel function) is not close to one of the nodes the com- 
parison between the solid line and the markers is good. 
However, from the inset we see that close to the node of 
Gi, we observe a relatively large discrepancy, which is 
obviously due to the correction from the 7r-periodic ef- 
fective potential given by Gi- The best parameter range 
to study the effect of the 7r-periodic effective potential is 
therefore to choose parameters such that G\ rs 0; i.e., for 
Jk (r) = when k is an integer. 

Figure 2 shows direct comparisons between numerical 
simulations (markers: a — 0.05 closed, a = 0.1 open) 
and our predictions (solid line), Arjk — 2| G 2 1 , for pa- 
rameter values leading to G\ = 0. The comparisons are 
performed at the smallest nonzero value of Y for which 
Jk{Y) — 0, and we show comparisons for k = (figure 
2a), k = 1 (figure 2b), k = 2 (figure 2c), and k = \ (fig- 
ure 2d) keeping Y constant for each figure. Note that in 
the latter case, k = |, Gi is always zero. We have here, 
arbitrarily, chosen Jo(Y = 2.4) ps 0.0. It is obvious that 
our comparisons demonstrate an excellent agreement be- 
tween simulations and prediction of the magnitude of the 
phase-locked region in n for all the different parameter 
values. The comparisons are performed in the frequency 
range between = 1 and £1 = 40 since driving frequen- 
cies smaller than 0,-2 typically lead to low stability of 
the phase locked states (the analysis is developed for high 
Q) and since phase-locking becomes impractical to iden- 
tify for frequencies larger than $1 = 30. Our data shows a 
slight trend of overestimating the locking range for large 
frequencies. We have identified this to be an artifact of 
numerically solving the pendulum equation with discrete 
time. Particularly, we have chosen the time step for the 
simulations to be a fraction of the period of the driving 
frequency, and thus, when looking for extremely small 
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phase-locked steps in the normalized IV characteristics, 
we are finding an artificial phase-locking of the dynamics 
to the temporal discretization. 

We have demonstrated that the 7r-periodic effective 
potentials can exist in the directly ac-driven pendulum 
and we have further given a quantitatively correct esti- 
mate of the significance of this effect. The magnitude of 
the 7r-periodic effective potential suggests that the phase- 
locking signature of the potential should be directly ob- 
servable not only in the driven pendulum jl5|,[l6), but 
also in ac-driven Josephson junctions. For relatively low 
driving frequencies, f2 ~ 2, we observe locking ranges 
in figure 2 of the order of Aij « 0.1, indicating that a 
standard dc current-voltage characteristic of a current 
ac-driven junction will indeed show a significant locking 
range where the usual Bessel function amplitude would 
suggest that locking is not possible. A particularly con- 
venient choice of parameters is to operate the system at 
the subharmonic resonance, k = 5, where G\ is always 
zero. 
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Appendix 

Let us consider a system of equations written in 
vector form 

X = f(x,T), 

where r = fit, x = (0,p), and / = [fx, f2). We will apply 
an averaging procedure as follows. 

Let x = x\ + Q,~ 1 hi{x\,T) be the first transformation 
with, yet undefined h\. This function is restricted to be 
periodic in r so that the new variables are close to the old 
ones uniformly in time. In the new variables the equation 
takes the form 

(I + Q^D^hx)^ + h lT = f(xx + Q^hi^r), r), 

where / is the unit matrix and D x denotes differentiation 
with respect to the elements in x. We obtain by Taylor 
expanding /, 

(I + 9r l D Xl hx)±x + hx T 
= f(xx,T) + n- 1 D x J(xx,T)h 1 +--- 

In order to eliminate the oscillatory part of f(xx, fit) we 
choose hi = {/(rzq, t)}_i ({</} = g - (g) and {g}-i is a 
mean-zero antiderivative of {<?}) to obtain 

±x = (/Xzij + n-^ifo.r.fi- 1 ), 

where R\ is polynomial in f} -1 . The second transfor- 
mation, given by x\ — x 2 + £l~ 2 /i2(:E2, t) with h 2 = 
{Rx(xx,t, 0)}_i, moves time-dependence to second order 

in rr 1 

±1 = (f)(x 2 ) + Q- 1 (R 1 )(x 2 ,0) + n- 2 R 2 (x 2 ,T,Q- 1 ). 

Continuing this procedure one can bring the system to 
the form 

Xn = (f){Xn) + n- 1 (R 1 ){ Xni 0) + 

■■■ + n- n+1 (Rn-x)(Xn,0) + n- n Rn(Xn,T, fT 1 ). 
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FIG. 1. (a) Simulated IV (?7, (</>)) characteristics of the ac-driven system described by Eq. ([j]). Parameters are: 57 = 3, 
a = 0.3, and r = 5. Relevant resonant k steps are indicated with arrows, and the range of phase-locking, Ar/k is indicated for 
the k = 1 step, (b) Magnitude, A771, of the phase-locked step as a function of e for parameters: Q — 3 and a — 0.1. Solid line 
represents the usual prediction of the Shapiro step, Ar/i = 2|Gi| for k = 1 and markers are results of numerical simulations. 
Inset shows details near the point Ji(F) = 0. 



4 




FIG. 2. Minimum of the locking range in r\ as a function of the driving frequency Q near the first node of Gi for nonzero T. 
Solid lines show the prediction, Arjk = 2|C?2 1 and markers are results of numerical simulations of Eq. ([[]). Open markers are for 
a = 0.1 and closed are for a = 0.05. (a) T « 2.4 and k = 0. (b) T » 3.8 and fc = 1. (c) T « 5.1 and fe = 2. (d)T« 2.4 and 
fc = I. 
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